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Abstract. - We study the effect of a Chern-Simons (CS) term in the phase structure of two 
different Abelian gauge theories. For the compact Maxwell-Chern-Simons theory, with the CS- 
term properly defined, we obtain that for values g = n/27r of the CS coupling with n = ±1, ±2, 
the theory is equivalent to a gas of closed loops with contact interaction, exhibiting a phase 
transition in the 3dXY universality class. We also employ Monte Carlo simulations to study 
the noncompact U{1) Abelian Ifiggs model with a CS term. Finite size scaling of the third 
moment of the action yields critical exponents a and v that vary continuously with the strength 
of the CS term, and a comparison with available analytical results is made. 



Gauge theories in 2 + 1 dimensions are frequently proposed as effective theories of strongly 
correlated electron systems in two spatial dimensions and zero temperature. Strong local 
constraints on the dynamics of lattice fermion systems are enforced by fluctuating compact 
gauge fields which exhibit topological defects in the form of space-time instantons [1]. It is 
conceivable that such effective gauge theories may feature phase transitions from a phase of 
bound instantons (deconfined phase) to a phase of unbound free instantons (confined phase). 
Such theories have been proposed as effective theories for high-temperature cuprate supercon- 
ductors, chiral spin liquids, and Mott insulators [2-4]. There is some hope that confinement- 
deconfincment transitions in 2 -t- 1 dimensions may shed light on quantum phenomena such as 
spin-charge separation and the breakdown of Fermi liquid theory in two spatial dimensions, 
and quantum phase transitions in Mott insulators. One central issue is how to describe ex- 
otic physical phenomena which do not comply with the Fermi-liquid paradigm such as the 
fractional quantum Hall effect, high- Tc superconductivity, and heavy fermion physics. 

In this paper we will consider two different (2 -I- l)-dimensional Abelian lattice gauge 
theories containing a Chern-Simons (CS) term [5]. The first of these theories is defined by 
a kind of lattice Maxwell-Chern-Simons Lagrangian. Polyakov has demonstrated that the 
compact Maxwell theory is permanently confined in 2 -I- 1 dimensions and does not exhibit 
any phase transition [1] . When this theory is coupled to bosonic matter fields with an integer 
charge g > 1, a deconfinement transition occurs [6,7]. The case where the bosonic matter 
fields are coupled to the fundamental charge is more controversial and at present it is not yet 
known if a deconfinement transition indeed occurs [8,9]. In this paper we will study the effect 
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of adding a CS-term to Polyakov theory. Such a theory was studied previously by a number 
of authors [10, 11]. In this paper we will use duality techniques to obtain new exact results on 
a lattice Abelian compact theory with a special type of CS-term. We will show that a phase 
transition does occur as a function of the gauge coupling for a fixed value of the CS-coupling, 
while confinement of electric charges is suppressed. 

The second theory considered in this paper is a non-compact lattice Abelian Higgs model 
with a CS-term. This theory has also been much studied in the past by many authors [12-15], 
often in connection with condensed matter systems. In this paper, we will compute the critical 
exponents of this theory as a function of the CS-coupling using Monte Carlo simulations. To 
our knowledge, this is the first time Monte Carlo simulations are employed on a lattice CS- 
theory. 

The lattice Maxwell-Chern-Simons theory with a compact gauge field used in this paper 
is defined by the partition function 



A = j^{A X A, - 2™,)' + f (A. - 27rN,) • (A x A, - 2™,), (1) 

where A^ is a periodic gauge field, and are integer fields taking care of the compactness 
of the theory and A is the lattice difference operator. The above Lagrangian is written in 
Villain-like form in order to allow the use of standard duality transformations. Note that it 
differs from the usual compact Maxwell-Chern-Simons Lagrangian [11], where the constraint 
A X Ni = rii is expected to hold. Physically, this constraint leads to a complete suppression 
of magnetic monopoles (which are instantons in 2 -I- 1 dimensions), since it implies A • rij = 0. 
Relaxation of this constraint incorporates the magnetic monopoles in the theory, since then 
A-iii = Qi, where Qi is an integer, just like in Polyakov's (2+l)-dimcnsional compact QED [1]. 
As we will see, removing the constraint A x N,; — iij implies new interesting physics. 

A Lagrangian similar to the one in Eq. |^ arises in effective descriptions of chiral spin 
states [16]. Due to compactness, Eq. Q is gauge invariant only if g = n/2TT, with n integer 
[10]. This is in contrast with the non-compact Abelian theory where g can be any real 
number. Furthermore, Eq. is also invariant under the integer gauge transformation Ai — *■ 
A., + 27rM,j, Ui ^ ni + Ax M^, and ^ -I- M^. 

Let us associate a vector field a with A x A — 27rn and b with A — 27rN by introducing 
the Lagrange multiplier vector fields A and a (in the following we omit sometimes lattice 
subindices to simplify the notation) 

C = + f a • b + a • (A X A - 2™ - a) + icr • (A - 27rN - b). (2) 

Next, we use the Poisson summation formula to replace A and cr by a new set of integer valued 
fields denoted L and S 

£ = ^g,^ + f a • b + iL ■ (A X A - a) + • (A - b). (3) 

Straightforward integration of A leads to the constraint A x L = S in the partition function. 
Summation over the field S followed by integration of the fields a and b yields 

£ = ^(AxL)2-^L.(AxL), (4) 

where we have written explicitly g — njlix. The crucial point to note is that when n — ±1, ±2, 
the integer CS-term in Eq, |0J does not contribute to the partition function. In such a case, 



E. SM0RGRAV, J. SmISETH, A. SUDB0 and F. S. NOGUEIRA: AbELIAN ChERN-SiMONS gauge THEORIES3 



the theory can be written in terms of a new integer field R = A x L to obtain 



where X]{r} denotes a constrained sum over closed vortex loops, i.e., with the constraint 
A • R = implied. The theory defined by Eq. ||SJ) is a theory of closed loops interacting 
through contact repulsion. This is the well known loop gas representation of the Abelian 
Higgs model with zero screening length which exhibits a loop-proliferation phase transition in 
the 3dXY universality class. There is a family of critical points given by w 0.33 x 87r^/n^, 
with the values n = ±1,±2. When n — ±2, the loop gas partition function ^ is equivalent 
to a non-compact abelian Higgs model in the so called "frozen" limit [17], which is the basis 
for the "inverted" 3dXY universality class in superconductors [18]. For n = ±1, on the 
other hand, Eq. ((SJ corresponds to a frozen superconductor where the charge of the Cooper 
pair / = 2e is fractionalized. Therefore, in this case the theory in ((SJ can be thought as 
corresponding to a "frozen" superconductor with charge //2 — e. 

The constraint A • R = implies that there are no monopoles in the spectrum of the theory 
with n — ±1, ±2. Therefore, there is no confinement of electric charges in the corresponding 
3dXY phase transition. This is entirely complementary to the theory without a CS-term: 
There, one has permanent confinement and no phase transition. 

In Ref. [11] a continuum version of a similar model is studied using a Hamiltonian approach 
via a variational analysis. There a ATy-like phase transition is also found, but only as a 
function of n. In contrast with our model, the compact Maxwell-Chern-Simons theory studied 
in Ref. [11] does not undergo any phase transition for fixed n. 

We next consider the Abelian Higgs model with a noncompact CS gauge field given in the 
Villain approximation by 



^(A0, - /A, - 2^n,)' + igA, • A x A, + A(A x A,)' 



(6) 



where A^ is a noncompact gauge field and 9i is a scalar phase field. Note that here the CS 
coupling g can be any real number, since now we are dealing with a non-compact Abelian 
gauge field. The topological defects of this model are vortex loops. An Abelian Higgs theory 
with a CS term has been proposed as an effective theory for the Laughlin state of fractional 
quantum Hall systems [3,12]. 

By introducing an auxiliary continuous field Vi, we can rewrite the action in Eq. as 



E 



-^v^ - ivi ■ {AO, - fAi - 27rni) + igA, ■ A x A, + A(A x A^)^ 



(7) 



Besides the usual gauge invariance, the action © is also invariant under the integer gauge 
transformation 9i 9i + 2TTli, rii + Ali, where k is an integer. This suggests that a gauge 

fixing on would allow us to extend the limit of integration for 6i over the whole real line. 
However, since we are dealing with integer fields, not all gauge fixings work consistently. The 
widely employed gauge fixing A • = does not work because it would lead to a decoupling 
of 6i from and consequently to a wrong loop-gas representation in the SdXY limit. Instead, 
we follow Ref. [19] and fix the axial gauge — 0. Within this gauge fixing we can integrate 
out 9i € {—oo, oo) in Eq. Q to obtain the constraint A ■ = 0. This constraint is solved by 
introducing a new vector field such that = A x h^. Next we perform a partial summation 
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in the term 27rmi • (A x h^) to give 2-Kih.i ■ (A x n^) and identify the integer field mi = A x rij 
as the vortex field. The end result is 



1 



— (A X h,;)-" - i/(A X h,;) • A, + 27ri • m, + igK, • A x A, + A(A x A,)" 
zp 



(8) 



where the constraint A • nii = is understood. Setting A = and integrating out A; yields 
a well known self-duality [15] at [3 ^ oo. 

After integrating out A; and and performing a lattice Fourier transform, we obtain 



-9' 



<i (¥'Q,Q-,+x)"+Q,Q-, 



(9) 



where QqQ-q = Q'^Q'^^, with — e**^/^ — e"*''^/^ being the Fourier representation of the 
symmetrized difference operator, and we have defined (p = X/g and x = /3/^/(2.9). V^'^ 
is the Fourier transform of the vortex-vortex interaction tensor. In the two limits g oo 
and f — ^ reduces to the interaction potential of vortices in the 3dXY model V^^^^ = 

i2TT^(3/QqQ.q)6f'-'. 

We next study the model Eq. @ using Monte Carlo (MC) simulations. Having real- 
valued vortex variables requires working in real-space. However, the action defined by Eq. 
is complex in real-space, 5* = Sr + iSi, leading to a complex transition-probability between 
various vortex-configurations in the Metropolis algorithm. This is analogous to the "sign 
problem" in quantum Monte Carlo simulations, which is often treated by using the absolute 
value of the probability p in the Metropolis update and sampling expectation values by (O) = 
{Os)^pJ {s)\p\, where s is the sign of the probability [20]. The expectation value of an operator 
O is in our case defined schematically by 

where denotes the complete set of eigenstates given by 5. Since the real part of the action 
Sr defines the same set of eigenstates, i.e. closed vortex loops, we may express {0)s in terms 
of expectation values {)sr 

where Zsfi is the partition function for the system defined by the real part of the action. The 
MC simulations are performed using e'^^ as the Boltzmann weight and expectation values are 
calculated from Eq. (fTT|l . 

The system size L and the lattice constant define the only length scales of the system. 
Critical properties are governed by the long range physics i.e. the q —^ limit. In this limit 
the real part of the potential has an effective screening length ^ x^/(2(/?x + !)■ When 
this screening length is of the order of L/2, the critical behavior will probably be a crossover 
to the 3dXY model with an infinite screening length. On the other hand, when AeS < a the 
critical behavior will experience a crossover towards a system with steric repulsion. Hence, 
MC simulations on finite lattices can only provide true critical exponents in a limited region 
of couphng space where a <C Aeg <C i/2. 
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Fig. 1- a.) {1 + a)/u from finite size scaling of A/3. The solid line denotes the well known 3dXY 
values and is included here for the sake of comparison, b) Critical exponent 1/ computed directly from 
the width of M3. c) a from combining M3 results for (1 + a)/!/ and l/i^. 



We have performed MC simulations on the model defined by 



(12) 



where the potential V^'^{r) is the inverse Fourier transform of Eq. A MC move is an 

attempt to insert a unitary closed vortex loop of random orientation, and the move is accepted 
or rejected according to the standard Metropolis algorithm using the real part of the action. 
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One sweep consists of traversing the lattice performing one MC move for each lattice site. 
The system with size L x L x L has periodic boundary conditions. Simulations with up to 10^ 
sweeps per coupling constant combined with Ferrenberg-Swendsen multihistogram reweighting 
are used to produce finite size scaling (FSS) plots. System sizes L g [4, 6, 8, 10, 12, 14, 16, 20, 24] 
have been used. 

To investigate the critical properties of the model, we compute critical exponents a and 
V. Recently we have proposed a FSS method based on the third moment of the action 
S [7], M3 {S — (S))^ >s, for which the asymptotically correct behavior is reached for 
accessibly small system sizes. The peak to peak value of this quantity scales with system size 
as and the width between the peaks scales as L~^/^ . In this way one can resolve a 

and V independently from one measurement without invoking hyperscaling [7] . 

The denominator \ {&''^')sr \ decreases with system size and coupling strength, contaminat- 
ing expectation values Eq. Hll() . Hence, as the strength of the CS term g increases, smaller 
system sizes are accessible for the FSS. Together with the limitations on the effective screen- 
ing length Aoff, this restricts the range in coupling space for which critical exponents can 
be calculated. We compute exponents along a line in coupling space corresponding to fixing 
/ = 1, A = I;, and tuning g so that Eq. (|ll|l is meaningful and crossover effects are absent. 
Along this line exponents can only be extracted from Ms-analysis for g < 0.5. We simulate 
system sizes L = 4, 6, 8, . . . where the maximum system size varies from i = 12 for g — 0.5 
to at least i = 24 for g = 0.01. The results of the FSS analysis are shown in Fig. We 
find non-universal exponents a and v approaching 3dXY values for large g and when g 0. 
Continuously varying critical exponents are a consequence of the vanishing of the renormaliza- 
tion group (RG) /3-function of the CS-coupling [13,21]. They are associated with a marginal 
operator, which in this case is just the CS-term. Critical exponents for the Abelian Higgs 
model with a CS-term have been obtained previously using RG methods [14,21]. Here we 
provide for the first time a MC calculation of these exponents. Due to the technical difficulties 
explained in Ref. [14], the RG calculations are more reliable for a generalized model with A^/2 
complex field components where N is large enough. Although our results are for N — 1, it 
is useful to compare at least qualitatively the results obtained here with those obtained from 
the 1/iV-expansion. In this case the critical exponent 1/ is given at order 1/A^ by [21] 
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ly = 1 — 



9 + ir \ ' ^'^^ 

where we have defined g = ^g/ir. The critical exponent v as given by Eq. 1)1 3|l exhibits 
a similar qualitative behavior in comparison with panel (b) in Fig. ^ having a maximum 
at some value of g. For smaller values of N , on the other hand, the RG treatment is in 
poor agreement with our numerical results, both from qualitative and quantitative points of 
view [14]. 

In summary, we have shown in two different gauge theories that the presence of a CS term 
changes dramatically their behavior. In both cases the analysis was entirely non-perturbative 
and based on duality arguments. 
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